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Abstract 

We investigate a strongly-correlated Bose gas in an optical lattice. Extending the standard- 
basis operator method developed by Haley and Erdos to a boson Hubbard model, we calculate 
excitation spectra in the superfluid phase, as well as in the Mott insulating phase, at T = 0. In 
the Mott phase, the excitation spectrum has a finite energy gap, reflecting the localized character 
of atoms. In the superfluid phase, the excitation spectrum is shown to have an itinerant-localized 
dual structure, where the gapless Bogoliubov mode (which describes the itinerant character of 
superfluid atoms) and a band with a finite energy gap coexist. We also show that the rf-tunneling 
current measurement would give a useful information about the duality of a strongly-correlated 
superfluid Bose gas near the superfluid-insulator transition. 

PACS numbers: 03.75.Lm, 03.75.Kk, 03.70.+k 
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Recently, a strongly-correlated Bose has been realized in an optical lattice P,0]. When the 
optical lattice potential is strong enough, the system can be described by a boson Hubbard 
model ^ 0] , where atoms are hopping between nearest-neighbor lattice sites, interacting with 
each other with the on-site repulsion U. In a cold atom gas, the hopping matrix element 
if: is tunable by adjusting the strength of laser light, so that one can realized the strongly- 
correlated regime characterized hy U/zt > 1 (where z is the number of nearest neighbor 
Sites). Theo.e«ca. wo.k on the boso HubbaM .odelg BBflBBynflhas 
predicted the superfluid-Mott insulator (SI) transition, which has been recently confirmed 
in Rb Bose gas The experimental result shows that the SI transition occurs at 

U/zt ^ 6 in a. cubic lattice, which is in good agreement with the theoretical prediction 
U/zt = 3 + 2\/2 = 5.830, O]. More recently, the optical lattice has been also applied to a 
cold Fermi gas^^, and the Brillouin zone has been observed. 

In the superfluid phase near the SI transition, while superfiuid atoms are itinerant, we 
can expect that the system would also have an insulating (localized) character. Indeed, an 
itinerant-localized dual character has been obtained in a fermion Hubbard model [14], where 
an itinerant heavy fermion band at the Fermi level and upper- and lower-Hubbard bands, 
describing the /oca/uec? character of fermions, coexist. At the metal-insulator transition, the 
itinerant heavy fermion band vanishes. In the Mott insulating phase, the upper and lower 
Hubbard bands only remains, that are separated by the so-called Hubbard gap Eg ~ U. 
Thus, it is a very interesting problem how this kind of duality can be seen in a strongly- 
correlated superfluid Bose gas. For this purpose, cold atom gases in optical lattices are very 
useful, because various physical parameters, such as the number density n of atoms and the 
interaction U/ zt, are tunable. Since the SI transition is absent in a weakly-interacting Bose 
gas, the duality is a characteristic phenomenon where the strong interaction between atoms 
plays an essential role. 

In this paper, we discuss excitation properties of a strongly-correlated Bose gas near the SI 
transition. Extending the standard basis operator method developed by Haley and Erdos^^ 
to the boson Hubbard model, we calculate the superfluid density of states at T = 0. This 
method has been recently used by Sheshadri[7], Kanobejiol and their coworkers to study the 
SI transition and the rf-spectrum in the Mott phase, respectively. In the superfluid phase 
near the SI transition, we show the coexistence of the itinerant gapless Bogoliubov mode 
and a band with a flnite threshold energy, which reflects the localized character of atomic 
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FIG. 1: Phase diagram of the boson Hubbard model. The SI phase boundary shown by the sohd 
(dashed) hne is given by (/i„) in ((111)) . The sohd circles at U/zt = 8 (a~h) are the positions 
where excitation spectrum and the dispersion relation are calculated in Figs. 2-4. fj,/zt at each 
position equals (a) ft- + (fi-^- — p,^)/3 (n = 1), (b) /i_ (n = 1), (c) 2/i_/3 (n = 0.89), (d) /i-/3 
(n = 0.74), (e) fx. + 2(/x+ - /i_)/3 (n = 1), (f) /x+ (n = 1), (g) /i+ + (C/ - A+)/3 (n = 1.16), and 
(h) /2+ + 2(C/ - /2+)/3 (n = 1.36). 

states. We also show that this novel phenomenon may be observed by the rf-tunneling 
measurement. 

We consider a superfluid Bose gas in a cubic optical lattice, described by the boson 
Hubbard model 



Here, q is a creation operator of a Bose atom at the i-th. site, and hi = c\ci. In the first 
term, t describes the nearest-neighbor hopping, where the summation («, j) is taken over 
the nearest-neighbor sites. U is the on-site repulsion, and /x is the chemical potential. For 
simplicity, we ignore effects of a trap potential in ((H). 

In the superfluid phase, we divide the Bose operator q into the condensate $ = (q) 
and the non-condensate 5ci = Ci — $. (We take $ to be real.) Substituting the expression 
Cj = $ + 6ci into ((H), we obtain H = ^ . hi — t j) SclScj, where hi = zt^"^ — zt^{ci + cj) + 
^hiifii — 1) — jJihi (where z = 6 is the number of nearest-neighbor sites of the cubic lattice). 
The BEC order parameter $ is determined by minimizing the ground state energy -Eq,=o of 
the on-site Hamiltonian h^^. Using the eigenstates \i.,a) of hi (with eigenenergies Ea)-, we 
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rewrite the Hamiltonian as 



a,al ,13,13' 



where T'^J^, ^ = {i, a\cl\i, a') {j, P\cj\j, P') . U^^, = \i,a){i,a'\ is referred to as the standard 
basis operators [l^, and ^L^^^' = ^a,a' ~ {^a,a')j where (■ ■ ■) is the expectation value in the 
ground state \i, a = 0). 

To calculate the excitation spectrum, it is convenient to introduce the single-particle 
Green's function, defined by gi,j{T) = — 'i6(r)([cj(r), cj]) [where 6(r) is the step function]. 
In the Lj^„,-representation, we obtain ^ij(r) = J2a,a',f3,i3'^'!i3',a,a'Ga,a',f3,t3'i^)^ where 

In the random phase approximation (RPA) in terms of the hopping t, the equation of motion 
for G is given by[3, Q|, in the energy {uj) and momentum (p) space, 

Sa,l3'Sa',f3Da,a' = —£pDa,a' ^g', 0,7,7' <J7,7', /3,/3' (P; ^) 

7.7' 

+ [U - {Ea' - Ea)]Ga,a',f3,P'{P,^^)- (4) 

Here, = —2t J2j=x y z '^'^^(Pj) describes the tunneling effect between sites (where the lattice 
constant is taken to be unity). D^a' = {La,a) - {La',a'), and f^',a,^^^' = T'J^a^n' + ^7,7',a',a' 
where j) is a nearest-neighbor pair. Substituting the solution of (0} into the single-particle 
Green's function (^(p, tu) in the energy and momentum space, we obtain 

I \ n(p,u;) 

1 - ep\l[p,uj) 



Here, 
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FIG. 2: Excitation spectrum N{uj) near the SI phase boundary at U/zt = 8. The value of fi in 
each panel is shown in Fig. 1. (a) and (e): Mott phase, (b) and (f): SI phase transition, (c) and 
(d): Superfluid phase with n < 1. (g) and (h): Superfluid phase with n > 1. We have introduced 
a small imaginary part T/zt = 0.001 to the energies, to smooth out the results. 

where = u + i6. We note that ya,i3 = {i, (y\ci\i, /?) is independent of the position i. 

In the Mott phase ($ = 0), the number density of atoms n per site is an integer, and the 



eigenstates \i, a) of hi are simply given by the number states, \i, n) 



c 



vacuum) with 



En = — + ^n{n — 1). Among them, |z, n = no) is the ground state when U{no — 1) < yU < 
Utiq. Noting that yn,n' = v^^n,n'-i, we obtain the Green's function in the Mott phase as 

Au{uj) 



1 - 6pAii{uj)^ 



where 
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n + 1 



n 
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Here, AEn+i,n = En+i — En = Un — fi is the excitation energy by adding an atom to the 
ground state when t = 0, while A_E„_i^„ = -En-i —En = fi — U{n — l) is the energy to extract 
an atom from the ground state. AEn-i.n may be regarded as the energy to add a hole to 
the ground state. The SI transition is determined by the condition when (jSI) has a gapless 
excitation at p = 0, which gives 3, Q, 0, 3| 

1 
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FIG. 3: Dispersion relation in the px direction, (b) SI phase boundary, (c) Superfluid phase. 



The lowest value of U/zt to realize the Mott phase equals U/zt = 3 + 2\/2 = 5.83 (n = 
3,0- We show the SI phase boundary given by ([111)) in Fig. 1. 
Now we consider excitation properties near the SI phase boundary. Figure 2 shows the 
density of states defined by 

N{uj) = --Y,lm[g{p,uj)]. (11) 
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In Fig. 2, the positive energy region describes particle-type excitations, and hole-type exci- 
tations are given by the spectrum in the negative energy region. We note that N{uj) < 
for uj < 0. [Since the Bose distribution function nsiuj) is negative for u < 0, the number of 
particles N{uj)nB{uj) at u is always positive.] In the Mott phase with n = 1 [panels (a) and 
(e)], a finite energy gap exists in particle and hole excitations. When t = 0, the energy gap 
Eg between particle and hole excitations equals Eg = AEiq + AEqi = U, so that Eg is found 
to correspond to the Hubbard gap in the fermion Hubbard model. Namely, the presence 
of this gap reflects the localized character of Bose atoms in the optical lattice. The finite 
widths of the particle and hole bands in Fig. 2 are due to the nearest neighbor hopping 
t. At the lower SI phase boundary [dashed line in Fig. 1], Fig. 2(b) shows that the hole 
excitations become gap less, which means that holes become itinerant. On the other hand, 
excitations in the particle branch {u > 0) still has a gap in Fig. 2(b). However, as one enters 
the superfluid phase [panels (c) and (d)], a gapless band suddenly appears in the particle 
branch. The resulting excitation spectrum in the particle branch has a dual structure, where 
the gapless excitation spectrum and the band with a finite threshold energy coexist. We 
note that this phenomenon can be also seen when one enters the superfluid phase through 
the upper SI phase boundary [solid line in Fig. 1]. In this shown in Figs. 2(e)-(h), 



particle excitations become gapless at the SI phase transition, and the dual band structure 
is obtained in the hole branch. 

Figure 3 shows the dispersion relation near the lower SI phase boundary. Although 
the spectral weights of the gapless particle and hole excitations are different in Fig. 2(c), 
their dispersion relations are symmetric with respect to u = 0. Their linear dispersions 
around Px = indicate that they are the Bogoliubov modes. At the SI transition, the 
gapless dispersion in the particle branch vanishes, and the hole dispersion becomes quadratic, 
reflecting the disappearance of superfluidity. 

The reason for the sudden appearance of the gapless spectrum in the superfiuid phase is 
that the macroscopic wave function of the superfiuid phase is a coherent state, involving the 
various number states \n) {n = 0, 1,2, ■ ■ ■). This leads to the hybridization of the particle 
and hole excitations in the superfiuid phase. As a result, the Bogoliubov excitation from the 
condensate is accompanied by both particle and hole excitations, which leads to the gapless 
particle and hole dispersions shown in Fig. 3 (dashed hues). 

We note that the gapless excitations in the superfiuid phase can be studied analytically 
in the strongly-correlated limit {U/zt ^ 1). In this extreme case, when n < 1, we may 
only take into account the two number states, \n = 0) and \n = 1), because of the strong 
repulsion. The resulting single-particle Green's function has the simple form 

^ (l-2nc)zf+(l-4rac)£p _ ^2ji _ 



n 



2 c<j_|_ — fip 

+ ^^E ^ 1 (12) 

2 u^ + Qp ' ^ ^ 

where Uc = $^ = n{l — n) is the condensate fraction. ()12j) shows the existence 
of the gapless dispersions in the particle and hole branches, given hj u = ±Qp = 



±a/ {zt + {2n — iyep){zt + Sp). In the long- wave length limit, we obtain u = ±'L>^p. Here, 
the sound velocity is given by v^/zt = $ a/2/3 = A^/2n^ /3jl7|. which is independent of the 
interaction U. This result is quite different from the usual Bogoliubov sound velocity in a 
weakly-interacting superfiuid Bose gas, given by = y^Un^jM oc VU (where M is the 
mass of an atom). We note that the spectral weight of the particle branch {u > 0) in p2|l 
vanishes at the SI phase boundary (ric = and = 1), as expected from Fig. 2(b). 

To observe excitation properties of a superfiuid Bose gas, the rf-tunneling current spec- 
troscopy is a powerful tooljlSf. In this measurement, atoms are transferred into an- 
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FIG. 4: Low energy rf-spectrum at T = 0. The contribution from the condensate /c at a; = is 
not shown in this figure. We have introduced a small imaginary part T / zt = 0.03 to the energies, 
to smooth out the results. 

other hyperfine state (described by a Bose operator 6p) by laser radiation. In the ro- 
tational wave approximation 3| ■, this process is described by the tunneling Hamiltonian 



7Ep 



e~'"^*&J,+q^Cp + h.c. 



where and ujl are the momentum and frequency 
of the laser light, respectively. 7 is the tunneling matrix element between the two hy- 
perfine states described by 6p and Cp. The Hamiltonian of the 6-state has the form 
Hb = X]p[^p + - l^bWpbp, where ^p = -2t ^^.^^ cosp^ + zt. oob and represent the 
threshold energy of the 6-band and the chemical potential, respectively. In this Hamiltonian, 
we have ignored the interaction in the 6-state and have assumed that the hopping matrix 
element t is the same as that of the c-state, for simplicity. In the first order perturbation in 
terms of Ht, we obtain the rf-tunneling current I {00) from the c-state to the 6-state as (we 
set = 0) 



(13) 



Here, u = ujl — Uh + fib is, a. renormalized detuning energy, and Ic = 27r7^?7,c describes the 
contribution from the condensate. The second term in involves information about the 
excitation spectrum shown in Fig. 2. Because of the step function in (fT!?|l . the second term 
is dominated by hole excitations at T = oQ|. 

Figure 4 shows the rf-spectrum near the Mott phase with n = 1. Near the lower SI 
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boundary [dashed line in Fig. 1], the excitation gap in the Mott phase simply vanishes at 
the SI transition [see panels (a)-(d)]. On the other hand, as shown in Figs. 4(e) and 4(f), the 
excitation gap in the Mott phase becomes larger near the upper SI phase boundary [solid 
line in Fig. 1]. The gap remains in the rf-spectrum even at the SI phase transition. In the 
superfluid phase, in addition to the band with a finite threshold energy, a gapless spectrum 
appears in the rf-spectrum [see Figs. 4(g) and 4(h)]. This result directly comes from the 
dual structure of the hole excitations shown in Figs. 2(g) and 2(h). Namely, the rf-spectrum 
in the superfluid phase with n < 1 is very different from that with n > 1. In particular, 
when n > 1, we can observe the dual band structure of a strongly-correlated superfluid Bose 
gas by the rf-tunneling measurement. 

To summarize, wc have discussed the excitation spectrum of a strongly-correlated Bose 
gas in an optical lattice. In the Mott phase, the Hubbard gap is obtained, reflecting the 
localized character of atoms. At the superfluid-Mott transition, particle or hole excitations 
become gapless while the other branch still has a flnite energy gap. However, the latter 
branch also becomes gapless in the superfluid phase due to the appearance of a gapless 
band. The resulting spectrum has a dual structure, where itinerant gapless excitations 
and the localized band with a finite threshold energy coexist. Wc have also shown that 
this dual structure of the excitation spectrum may be observed by the rf-tunneling current 
measurement. Since the duality is absent in a weakly-interacting Bose gas, the observation of 
this phenomenon would be very useful in understanding the physics of a strongly-correlated 
superfluid Bose system. 

This work was supported by a Grant-in- Aid for Scientiflc research from the Ministry of 
Education of Japan. 



[1] M. Greiner, et al., Nature 415, 39 (2002). 

[2] T. St5fcrlc, ct al., Phys. Rev. Lett. 92, 130403 (2004). 
[3] M. Fisher, et al, Phys. Rev. B 40, 546 (1989). 

[4] S. Sachdev, Quantum Phase Transitions (Cambridge University Press, Cambridge, 1999) 
Chap. 10. 

[5] D. Rokhsar, et al., Phys. Rev. B 44, 10328 (1991). 



9 



3] W. Krauth, et al., Phys. Rev. Lett. 67, 2307 (1991). 

7] K. Sheshadri, et al., Europhys. Lett. 22, 257 (1993). 

3] J. Freericks, et al., Europhys. Lett. 26, 545 (1994). 

3] D. van Oosten, et al., Phys. Rev. A 63, 053601 (2001). 

3] D. Dickerscheid, et al., Phys. Rev. A 68, 043623 (2003). 

1] R. Roth and K. Burnett, Phys. Rev. A 67, 031602 (2003). 

2] K. Sengupta, et al., Phys. Rev. A 71, 033629 (2005). 

3] M. Kohl, et al., Phys. Rev. Lett. 94, 080403 (2005). 

4] See, for exmaple, R. Bulla, Phys. Rev. Lett. 83, 136 (1999). 

5] S. Haley, and P. Erdos, Phys. Rev. B 5, 1106 (1972). 

3] S. Kanobe, et al., |cond-mat/ 0407229| 

7] In a d-dimensional lattice, we obtain v^p/zt = which agrees with Ref.[3| for a two- 

dimensional square lattice, giving V/p/zt = <I>. 

[18] D. Luxat, and A. Griffin, Phys. Rev. A 65, 043618 (2002). 



10 



